The superconducting phase of Calcium under the pressure at 200 GPa: the 

strong-coupling description 
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The thermodynamic parameters of the superconducting state in Calcium under the pressure 
at 200 GPa have been determined. The numerical analysis by using the Eliashberg equations in 
the mixed representation has been conducted. It has been stated, that the critical temperature 
(Tc) decreases from 36.15 K to 20.79 K dependently on the assumed value of the Coulomb 
pseudopotential (p* £ (0.1,0.3)). Next, the order parameter near the temperature of zero Kelvin 
(A (0)) has been obtained. It has been proven, that the dimensionless ratio 2A (0) jksTc decreases 
from 4.25 to 3.90 together with the growth of \x* . Finally, the ratio of the electron effective mass 
to the electron bare mass (m*/m e ) has been calculated. It has been shown, that m%/m e takes the 
high value in the whole range of the superconducting phase's existence, and its maximum is equal 
to 2.23 for T = T c . 
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I. INTRODUCTION 

In Calcium under the influence of the high pressure (p), 
the whole sequence of the structural phase transitions is 
induced. 

In the normal conditions Calcium crystallizes in the fee 
structure (Ca-I). For the pressure at 20 GPa, the tran- 
sition to the structure bec (Ca-II) has been observed. 
The next change appears for the pressure's value at 32 
GPa, above which the structure sc (Ca-III) is realized 
In 2005 Yabuuchi et al. have found another two crystal 
phases: Ca-IV and Ca-V 0, 0. It has been stated, that 
the phase Ca-IV is stable from 119 GPa; above the pres- 
sure of 143 GPa appears Ca-V. The phases Ca-IV and 
Ca-V have been indentified experimentally by Fujihisa et 
al. [!]• The following assignment has been proposed: Ca- 
IV with the P4i2i2 structure and Ca-V with the Cmca 
structure. In the range of pressures from 158 GPa to 
207 GPa, Nakamoto et al. have discovered the phase Ca- 
VI (the Pnma structure) Above, Sakata et al. have 
reported existence of the host-guest structure (Ca-VII) 
1- 

The results achieved by using the ab initio methods 
only partially agree with the experimental scheme of the 
structural phase transitions. In particular, Ishikawa et al. 
and Arapan et al. have reproduced the appearance of the 
low pressure structures (fee and bec) 0], [1]. However, 
the values of the pressure, at which the structural transi- 
tions appear, differ from the experimental values. Above 
the pressure at 36 GPa, the sequences of the succes- 
sive structural transitions, explicitly differ between them- 
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selves and do not reproduce the experimental scheme. In 
the context of the essential theoretical results, we have 
mentioned the paper of Yao et al. Q . The authors were 
the ones who identified the phase Ca-IV as the structure 
Pnma and the phase Ca-V as Cmca. 

In the range of the very high pressures (from 135 GPa 
to 495 GPa), the calculations of Ishikawa et al. suggest 
the stability of the structure I4:/mcm(00j). However, the 
similar value of the enthalpy has the Pnma structure (at 
last for the lower pressures from the discussed range) [7[ . 
It has to be underlined, that the stability of the Pnma 
structure from 158 GPa to 180 GPa is suggested by the 
paper of Aftabuzzaman and Islam Also, the results 
presented by Yin et al. prove especially low values of the 
Pnma enthalpy (p £ (135, 220) GPa) [n|. 

Besides the very complex sequence of the structural 
transitions, which to the present day has not been finally 
determined, Calcium is being characterized with the in- 
teresting superconducting properties. In particular, the 
very high values of the critical temperature among the 
simple metals. 

The possibility of the superconducting state's exis- 
tence in Calcium has been reported in 1981 by Dunn 
and Bundy [l^ . In 1996, the definite experiment has 
been performed by Okada et al. [l3|]. After ten years, 
the studies of Okada have been repeated by Yabuuchi 
et al. [14j. The researchers determined the dependence 
of the critical temperature on the pressure, at the same 
time ascertaining, that Tc appreciably exceeds the values 
determined by Okada. 

According to Yabuuchi et al., the superconducting 
state in Calcium is induced in the following crystal 
phases: Ca-III. Ca-IV and Ca-V [l4j|. In particular, in the 
phase Ca-III the critical temperature quickly increases 
together with the pressure's growth (from about 3 K for 
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p = 58 GPa to 23 K for p = 113 GPa). In the phase 
Ca-IV, the critical temperature is the subject of the sat- 
uration, while in the phase Ca-V it takes the value of 
25 K at 161 GPa. In 2011, Sakata et al. have reported 
T c = 29 K for p = 216 GPa @; the highest observed 
value of the critical temperature among all elements. 

The above experimental results inspired us to estimate 
the basic thermodynamic parameters of the Ca supercon- 
ducting state in the range of the very high pressures. In 
particular, we assume: p — 200 GPa. The calculations 
will be conducted in the framework of the Eliashbcrg 
formalism, with the aid of the Eliashberg function calcu- 
lated by Yin et al. in the paper [Tl| (Pnma structure). 
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II. THE ELIASHBERG EQUATIONS 



The Eliashberg equations have been derived in order to 
quantitative description of the superconducting state for 
the intermediate and strong coupling between electrons 
and phonons. 

In appliance to the approach used in the BCS theory 
[HI , the Eliashberg equations enable the detailed consid- 
eration of the complicated form of the electron-phonon 
interaction. In the result, the Eliashberg set allows ex- 
actly to estimate the value of the critical temperature, 
the order parameter and the electron effective mass, if 
the value of the Coulomb pscudopotential is known. 

The Eliashberg equations can be written on the imagi- 
nary axis, on the real axis or in the mixed representation 
[l6l | . From the numerical point of view, the most conve- 
nient in application is the mixed representation, because 
the rather simple algorithms can be used [l?}, (l8j . 
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FIG. 1: The dependence of the criticai temperature on the 
Couiomb pseudopotential. In the inset, we have shown the 
ratios r u and r\ as the function of the Coulomb pseudopo- 
tential. 
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where: K , >— 

The symbols: (j)(co), Z (co), (<p n = (f>(iuj n ) and Z n = 
Z (ico n )) represent the order parameter function and the 
wave function renormalization factor on the real (imag- 
inary) axis, respectively; co n is the n-th Matsubara fre- 
quency: co n = (tt/P) (2n — 1). The inverse temperature 
is given by: fi = (fc^T) -1 , and ks is the Boltzmann 
constant. The order parameter is defined as: A = cj>/Z. 

In the framework of the Eliashberg formalism, 
the pairing kernel for the electron-phonon interac- 
tion is determined in the following way: A (z) = 

2 / o nmax dQ,-^^a 2 F (O). The symbol a 2 F (fi) denotes 
the Eliashberg function (Tlj : the value of the maximum 
phonon frequency (ri max ) is equal to 78.11 meV. 

The depairing Coulomb interaction is modeled para- 
metrically with the aid of the Coulomb pseudopotential 
/x* . Symbol 9 denotes the Heaviside unit function and co c 
is the cut-off frequency (co c = 30 max ). The Bose-Einstein 
and Fcrmi-Dirac functions are represented by N (uj) and 
/ (to), respectively. 

The Eliashberg equations have been solved numerically 
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FIG. 2: The real and imaginary part of the order parameter 
on the real axis for selected values of the temperature. The 
rescaled Eliashberg function has been also plotted. 
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FIG. 3: The form of the order parameter on the complex plane 
for selected values of the temperature. The results obtained 
for uj £ (0, fimax) by using the lines with empty symbols have 
been denoted; the solid lines represent results achieved for 

W G (f2max,W c ). 



for 2201 Matsubara frequencies {A I = 1100). The stabil- 
ity of the solutions has been achieved for the tempera- 
tures greater or equal to Tq = 8.12 K. Additional infor- 
mation on the means of the analysis of the Eliashberg 
equations, the reader can find in the papers: [2(|, plj . 
and (H. 



III. THE RESULTS 

In the Eliashberg formalism, the critical temperature 
depends on the form of the Eliashberg function and the 
value of the Coulomb pscudopotcntial. In the opposition 
to the Eliashberg function, the parameter /i* is really 



hard to calculate using the ab initio methods. According 
to the above, we have taken into consideration the very 
wide range of the Coulomb pseudopotential's values e.g. 
fx* € (0.1,0.3). The significant suggestion, that not only 
low values of fx* should be considered, are the results 
obtained for Calcium under pressure at 120 GPa, where 
fx* = 0.215 has been noted [22| . 

In Fig. [TJ we have presented the critical temperature as 
the function of /x*. The exact results obtained by using 
the Eliashberg equations are represented by the filled cir- 
cles. It can be easily noticed, that together with increase 
of the parameter fx*, the critical temperature decreases 
from the value of 36.15 K to 20.79 K. The above result 
proves, that even for the very large value of Coulomb 
pscudopotcntial, Tc achieves high values. We notice, 
that the range of the calculated values of Tc agrees qual- 
itatively with the recently experimental data. In par- 
ticular, Sakata et al. have shown that: Tc = 29 K for 
p = 216 GPa 0. 

From the mathematical point of view the determina- 
tion of the critical temperature with the aid of the Eliash- 
berg equations is the complicated and time consuming 
issue. For that reason, we have presented the analytical 
formula, which exactly reproduces the numerical results 
in opposition to the classical expressions of Allen-Dynes 
or McMillan (see Fig.[TJ HU, In order to do that, 

the Allen-Dynes formula with the newly fitted parame- 
ters has been used. We notice, that the fitted parameters 
have been determined on the basis of the 250 exact val- 
ues of Tc (a**), and the least squares method. The result 
takes the form: 



ksTc = /i/2Y^ ex P 



-1.03(1 + A) 



A-/i*(l + 0.06A) 



(3) 



where the strong-coupling correction function (/i) and 
the shape correction function (/ 2 ) are given by the 



expressions: 



/l = 



1 
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The symbol A denotes the clcctron-phonon 



coupling constant: A = 2 J ' max d£l a ■ The quan- 
tity cj2 represents the second moment of the normal- 
ized weight function: u>2 = f fo**""' d£la 2 F (Q.) fl, and 
win is called the logarithmic phonon frequency: ui\ n = 

In particular, for Cal- 



cxp 



_|j?-dn^]n(n)j 

cium under the pressure at 200 GPa, the following values 
have been obtained: A = 1.23, -^/wi = 29.98 meV, and 
uj\ n = 35.92 meV. The functions Ai and A2 are defined as: 

Ai = 2.6 (1 + T8/1*), and A 2 = 0.092 (1 - 150//*) . 

The modified Allen-Dynes (mAD) formula possesses 
the fitted parameters, which are sharply different than 
the parameters included in the classical expression de- 
rived by Allen and Dynes (AD). The above result is con- 
nected with the fact, that the effective phonon frequency 
Qwefl^AD) an( I ^ nc effective coupling constant ([A c ff] AD ) in 
classical Allen-Dynes formula well reconstruct the exact 
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value of Tc only for low values of the Coulomb pseu- 

dopotential (fi* — > 0.1). We notice, that the parameters 

_ i 

Wed and A e ff have been defined as: ksTc = 0J e se Acff ■ 
In the case of high values of /i*, the classical parame- 
terization overestimates the effective phonon frequency 
and much underestimates the effective coupling con- 
stant. The inset in the Fig. 1 presents the detailed 
dependences of the ratios r u = [w e ff] mAD / [weff]AD an< ^ 

r\ = [Acff] mA D / [ A cff] A D on M*- 

Finally, we notice that the classical expressions for Tc 
should be applied with the great carefulness, because 
even for low-temperature superconductors like Al or Pb, 
the difference between analytical and numerical Eliash- 
berg results are noticeable |25j . 

In Fig. [21 we have presented the order parameter on 
the real axis for the selected values of the temperature 
and /x* = 0.1. We see, that for the low frequencies, the 
non-zero value is taken only by the real part of A(w). 
The obtained result proves, that in the considered range 
of frequencies the damping effects not exist [2(| . For the 
higher values of the frequency (ui ~ 30 meV) one can ob- 
serve the maximum of Rc[A (u>)], which is induced by the 
area of the characteristic peaks in the Eliashbcrg func- 
tion. From the physical point of view, the above fact in- 
dicates that the relevant growth of Re[A (a;)] exists in the 
frequencies' area in which the electron-phonon coupling is 
exceptionally strong. Additionally, it should be marked 
out, that in the considered case the function Im[A (uj)] in- 
creases monotonically together with the growth of uj. For 
the frequencies higher than ~ 30 meV, the real part of 
the order parameter decreases together with the growth 
of ui due to the drop of the Eliashbcrg function's value. 
In the last step let us turn our attention to the fact, that 
the growth of the Coulomb pseudopotential causes only 
the significant decrease of the order parameter while the 
shapes of Re [A (uj)] and Im[A (uj)] remain very similar. 

In Fig.[3J we have plotted the order parameter on the 
complex plane. It can be easily noticed, that the values 
of A (uj) make the characteristic spirals with the radius 
that decreases together with the growth of the tempera- 
ture. The obtained result allows to characterize the ef- 
fective electron-electron interaction. In particular, for 
Re[A (uj)]> the effective interaction is pairing, whereas 
for Re[A(w)]< the depairing processes are in major. 
In the case when /i* = 0.1, the area of the frequencies in 
which the electron-electron interaction is attractive cov- 
ers the range of the frequencies for which the Eliashbcrg 
function is definite. As the result of the /i* increasing, 
the frequencies that correspond to the positive value of 
Rc[A (a;)] become narrow from the side of fi m a X ; however 
such effect is not exceptionally strong. 

On the basis of the obtained results, we have calcu- 
lated the physical value of the order parameter for the 
given temperature (A(T)). In order to achieve that, 
one should use the following equation [l6| : A (T) = 
Re [A(uj = A (T))]. From the physical point of view, 
the most interesting is the value of the order parameter 
close to the temperature of zero Kelvin (A (0) ~ A (To)), 
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FIG. 4: The dependence of the ratio Ri on the value of the 
Coulomb pseudopotential. In the inset, the influence of n* on 
the value of the parameter A (0) has been presented. 
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FIG. 5: The form of the wave function renormalization factor 
on the complex plane for selected values of the temperature 
and fi* = 0.1. The results for uj £ (0, f2 max ) by using the 
lines with empty symbols have been denoted; the solid lines 
represent the results for uj £ (f2 ma x,kJ c ). In the inset, we have 
presented the dependence of the ratio m*/m e on the temper- 
ature for selected values of the Coulomb pseudopotential. 



because it allows to determine the characteristic ratio: 
i?i = 2A (0) /ksTc- Let us notice, that in the framework 
of the BCS theory, the parameter Ri takes the universal 
value equal to 3.53 [l5l ]. 

In Fig. HI we have presented the results achieved for 
Calcium, where the dependence of R\ on the Coulomb 
pseudopotential has been plotted. We see, that the 
parameter R\ is always bigger than [-Ri] BCS and only 
slightly decreases together with the increase of fi* (from 
4.25 to 3.90). We notice, that the weak dependence of 
the ratio i?i on /.i* is connected with the fact that the 
Coulomb pseudopotential in the comparable way lowers 
the value of A (0) and Tc- In particular, the open de- 
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pendence of A (0) on fx* in the Fig's. 4 inset have been 
plotted, whereas the shape of the function Tq (fJ.*) in Fig. 
Q]has been shown. 

From the physical point of view, the high values 
of the ratio [i?i] Ca are the consequence of the exis- 
tence of the strong-coupling and retardation effects, 
which are omitted by the simple BCS theory. In 
the framework of the Eliashbcrg approach, these ef- 
fects can be characterized with the aid of the param- 
eter ksTc/^xu- In the weak-coupling limit, one can 
assume: [k B T c /uJi n ] BCS = 0. In the case of Cal- 
cium, we have obtained: [ksTc /u)\ n ] u*=o i — 0.104, and 
[ksTc / uj\ n ] ^ =0 3 — 0.059. Above results clearly differ 
from the BCS approximation. 

Finally, we notice that the exact value of /z* for Cal- 
cium under the pressure at 200 GPa is unknown. How- 
ever, on the basis of the modified Allen-Dynes formula 
and the experimental result presented recently by Sakata 
et al. (T c = 29 K for p = 216 GPa) 0, one can obtain 
[M*]p= 2 i6GPa - 0- 177 - I f wc assume that: [M*] p=200GPa ~ 

[M*] p= 216GPa> thG ValuC 01 [ i? l]p=200GPa CC l Uals ab ° ut 4 -°' 

In the Eliashbcrg formalism the dependence of the elec- 
tron effective mass (m*) on the temperature is deter- 
mined by the wave function renormalization factor. In 
particular, we have: m* = Re [Z (0)] m e , where m e de- 
notes the electron bare mass. 

In Fig. [51 we have plotted the values of Z (lj) on the 
complex plane in order to investigate the influence of the 
temperature on the form of the wave function renormal- 
ization factor. On the basis of the obtained results, it 
has been stated, that the wave function renormalization 
factor weakly depends on the temperature in comparison 
with the order parameter. We notice, that the increasing 
value of the Coulomb pseudopotential does not change 
significantly the course of Z (ui). 

In the Fig's. 5 inset, the dependence of the ratio 
m*/m e on the temperature for the selected values of the 
Coulomb pseudopotential has been shown. The achieved 
results prove, that the electron effective mass takes the 
very high value in the whole range of the superconducting 
phase's existence, and m* has the maximum for T = Tc, 
where m* e /m e = 2.23. 

IV. SUMMARY 

In the paper, we have obtained the basic thermody- 
namic parameters of the superconducting state in Cal- 



cium under the pressure at 200 GPa. It has been stated, 
that dependently on the assumed value of the Coulomb 
pseudopotential (/i* G (0.1,0.3)), the critical tempera- 
ture decreases in the range from 36.15 K to 20.79 K. 
This result is in agreement with the recently experimen- 
tal data obtained by Sakata et al. ( Tc = 29 K for 
p = 216 GPa) 0. 

Next, the dependence of the order parameter A (0) on 
H* has been determined. The results allowed to deter- 
mine the value of the ratio R\. We have affirmed, that 
the parameter R\ is essentially larger than [i?i] B cg, and 
slightly decreases with the growth of the Coulomb pseu- 
dopotential from 4.25 to 3.90. 

In the last step, the electron effective mass has been 
determined. We have shown, that m* is large and takes 
its maximum for T — Tc- In particular: [w*] max = 
2.23m e . 

In the paper, the detailed characteristic of the order 
parameter on the real axis has been also presented. It has 
been proven, that the courses of Re[A (u>)] and Im[A (lj)] 
are clearly correlated with the shape of the Eliashbcrg 
function. Additionally, the plot of the order parameter 
on the complex plane enabled the conclusion, that for 
/j* =0.1, the effective electron-electron interaction is at- 
tractive in the frequencies' range, where the Eliashbcrg 
function is definite. The increase of /i* causes only the 
insignificant narrowing of the frequencies' range from the 
side of S! max - 
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